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We discuss the imphcations of assuming a four-zero Yukawa texture for the properties of the 
charged Higgs boson within the context of the general 2-Higgs Doublet Model of Type III. 
We begin by presenting a detailed analysis of the charged Higgs boson couplings with heavy 
quarks and the resulting pattern for its decays. The production of charged Higgs bosons 
is also sensitive to the modifications of its couplings, so that we also evaluate the resulting 
effects on the top decay t hH^ as well as on 'direct' cb + c. c. and 'indirect' 

qq, gg — s- thH^ + c. c. production. Significant scope exists at the Large Hadron Collider for 
several production and decay channels combined to enable one to distinguish between 
such a model and alternative 2-Higgs doublet scenarios. 
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I. INTRODUCTION 

Detecting a charged Higgs boson during the imminent Large Hadron CoUider (LHC) experi- 
mental running would constitute a clear evidence of physics beyond the Standard Model (SM) 



Charged Higgs bosons appear in many wel 
has been widely studied over the years 



motivated extensions of the SM, whose phenomenology 
ali. In particular, 2-Higgs Doublet Models (2HDMs), 



IS 



can be considered as a prototype of 



in both Supersymmetry (SUSY) and non-SUSY versions 
a Higgs sector that includes a charged Higgs boson (H^). It is expected that the LHC will allow us 
to test the mechanism of Electro- Weak Symmetry Breaking (EWSB) and, in particular, to probe 



y-interacting 
and general 
7|, or whether 



the properties of charged Higgs bosons, which represent a unique probe of a weak! 
theory, as is the case of the Minimal Supersymmetric Standard Model (MSSM) Isl 
2HDMs of Type I, II, HI and IV (2HDM-I, 2HDM-I, 2HDM-III and 2HDM-IV) 
strongly-interacting scenarios are instead realized, like in the old Technicolor models or similar ones 
discussed more recently Q]. Ultimately, while many analyses in this direction can be carried out 
at the LHC, it will be a future International Linear Collider (ILC) 191 or Compact Linear Collider 

n 

(CLIC) which will have the definite word about exactly which mechanism of mass generation 
and which realization of it occurs in Nature. 

The 2HDM-II has been quite attractive to date, in part because it coincides with the Higgs 
sector of the MSSM, wherein each Higgs doublet couples to the u- or d-type fermions separately^. 
However, this is only valid at tree-level jjj]. When radiative effects are included, it turns out 
that the MSSM Higgs sector corresponds to the most general version of the 2HDM, namely the 
2HDM-III, whereby both Higgs fields couple to both quarks and leptons. Thus, we can consider 
the 2HDM-III as a generic description of physics at a higher scale (of order TeV or maybe even 
higher), whose low energy imprints are reflected in the Yukawa coupling structure. With this idea in 
mind, some of us have presented a detailed study of the 2HDM-III Yukawa Lagrangian [l^, under 
the assumption of a specific texture pattern [l4], which generalizes the original model of Ref. [l^. 
Phenomenological implications of this model for the neutral Higgs sector, including Lepton Flavour 
Violation (LFV) and /or Flavour Changing Neutral Currents (FCNCs) have been presented in a 



previous work [16(]. Here, we are interested in extending such an approach to investigate charged 
Higgs boson phenomenology: namely, we want to study the implications of this Yukawa texture 
for the charged Higgs boson properties (masses and couplings) and discuss in detail the resulting 



^ Notice that there exist significant differences between the 2HDM-II and MSSM though, when it comes to their 
mass/coupling configurations and possible Higgs signals pH ]. 
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pattern of charged Higgs boson decays and main production reactions at the LHC. 

Decays of charged Higgs bosons have been studied in the hterature, including the radiative 

, n n 

modes W Z [17], mostly within the context of the 2HDM-II or its SUSY incarnation 



(i.e., the MSSM), but also by using an effective Lagrangian extension of the 2HDM 



181] and. 



more recently, within an extension of the MSSM with one Complex Higgs Triplet (MSSM+ICHT) 



211 ] and, more 
21. 



Charged Higgs boson production at hadron colliders was studied long ago 
recently, systematic calculations of production processes at the LHC have been presented 

Current bounds on the mass of a charged Higgs boson have been obtained at Tevatron, by 
studying the top decay t — > bH~^, which already eliminates large regions of the parameter space 



231 ]. whereas LEP2 bounds imply that, approximately, > 80 GeV 



24 



251 ]. rather model 



independently. Concerning theoretical limits, tree-level unitarity bounds on the 2HDM Higgs 
masses have been studied in generic 2HDMs and in particular an upper limit for the charged Higgs 



26]. 



mass of 800 GeV or so can be obtained, according to the results of Ref. 

This paper is organized as follows. In section II, we discuss the Higgs- Yukawa sector of the 
2HDM-III, in particular, we derive the expressions for the charged Higgs boson couplings to heavy 
fermions. Then, in section HI, we derive the expressions for the decays H'^ fifj and numerical 
results are presented for some 2HDM-III scenarios, defined for phenomenological purposes. A 
discussion of the main production mechanisms at the LHC is presented in section IV. These include 
the top decay t — > bH^ as well as s-channel production of charged Higgs bosons through cb{ch)- 
fusion [2^ and the multi-body more qq, gg — > tbH~+ c.c. (charge conjugated). These mechanisms 
depend crucially on the parameters of the underlying model and large deviations should be expected 
in the 2HDM-III with respect to the 2HDM-II. Actual LHC event rates are given in section V. 
Finally, we summarize our results and present the conclusions in section VI. Notice that in carrying 
out this plan, unlike other references 



2S 



2g|, where the 2HDM-II and the 2HDM-III appear as 
different structures, we shall consider here that, under certain limits, the 2HDM-III reduces to the 
2HDM-II and, therefore, that the properties of the charged Higgs bosons change continuously from 
one model to the other. 



II. THE CHARGED HIGGS BOSON LAGRANGIAN AND THE FERMIONIC 

COUPLINGS 

We shall follow Refs. [l^ . where a specific four-zero texture has been implemented for the 
Yukawa matrices within the 2HDM-III. This allows one to express the couplings of the neutral and 
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charged Higgs bosons in terms of the fermion masses, Cabibbo-Kobayashi-Maskawa (CKM) mixing 
angles and certain dimensionless parameters, which are to be bounded by current experimental 
constraints. Thus, in order to derive the interactions of the charged Higgs boson, the Yukawa 
Lagrangian is written as follows: 



Cy = Y^'Ql^^ur + Y^Ql^2Ur + YtQL^^dR + Y^QMr, 



(1) 



where $1^2 = ('^i'2' '?^i,2)"^ refer to the two Higgs doublets, $i_2 = io"2^i,2) Ql denotes the left- 
handed fermion doublet, ur and cLr are the right-handed fermions singlets and, finally, Y^2 denote 
the (3 X 3) Yukawa matrices. Similarly, one can write the corresponding Lagrangian for leptons. 

After spontaneous EWSB and including the diagonalizing matrices for quarks and Higgs 
bosons^, the interactions of the charge Higgs boson with quark pairs acquire the following 
form: 



H+ 



2V2MW 




Ui { {VcKM)ii tan j3 nidi Sij - sec (3 



+ 



cot P m„ . 6ii - CSC P 



CKMjjZ 



tan P rudi ~ 



cot P ruu, Sii - CSC P 




(2) 



v ^ 

where Vckm denotes the mixing matrices of the quark sector (and similarly for the leptons). The 
term proportional to 6ij corresponds to the contribution that would arise within the 2HDM-n, 
while the terms proportional to Y2 and Y^ denote the new contributions from the 2HDM-ni. 



These contributions, depend on the rotated matrices: Yn = OT Pq Yn Pq Oq {n = 1 when q 



and n = 2 when q = d ), where Oq is the diagonalizing matrix, while Pq includes the phases of 
the Yukawa matrix. In order to evaluate Yn we shall consider that all Yukawa matrices have the 
Hermitian four-zero texture form |il4|], and the quark masses have the same form, which are given 
by: 



Ml 



^ Cq \ 

C* Bq Bq 



Bl A. 



{q = u,d). 



(3) 



1/ 



The details of both diagonahzations are presented in Ref. [IS 
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This is called a four-zero texture because one assumes that the Yukawa matrices are Hermitian, 
therefore each u and d type Yukawa matrix contains two independent zeros. According to current 
analyzes this type of texture satisfies the experimental constraints and at the same time it permits 
to derive analytical expressions for the Higgs boson fermion couplings. 

To diagonalize M"^, we use the matrices Og and Pg, in the following way [l^ : 



M'^ =0'^PgMlP^Og. 



(4) 



Then, one can derive a better approximation for the product OT Pg Yn Pg Og, expressing the rotated 



matrix Yn, in the form 



V 



(5) 



In order to perform our phenomenological study, we find it convenient to rewrite the Lagrangian 
given in Eq. ([2]) in terms of the coefficients [xn]jj, as follows: 

3 



^ ^ sec/3 ^ J 

tan I3mdi dij ■j=- y/mj/mjj Xij 



+ 



cot f3 rriu^ - ^m„,m„, Xu 



V2 



V2 

(VcKM)zj 



CKMJil 



^ ^ sec/3 ^ J 

tan 13 nidj dij ^ ^md,md^ Xij 



V2 



7 



cot l3 ruu^ 5ii - '^^^ ^Jmu^mu^ Xu 



{VcKM)ij r } dj H 



(6) 



where we have redefined [Xi]ij = X^j aiid [xfljj — xfj- Then, from Eq. the couplings UidjH' 
and UidjH^ are given by: 



(7) 



where Sij and Pij are defined as: 



Sij — y^(VcKM)» 



1=1 

+ 

3 



^ ^ sec/3 , J 

tan [3 rudi ^ij 7^ ^/mJjnJJXij 



V2 



^ r, r esc/? „■ 

cot 13 niu, Oil — y/mu.mui Xu 



V2 



(VcKM)/j, 



= 



CKM)U 



1=1 



^ J. sec/3 , J 

tan 13 nidi ^ij 7^ ^/rrW^ Xij 



V2 



^ r, r csc/3 „■ 

cot (3 rriu, 6u - —j^ ^Jmu^niui Xu 



(VcKM)/j- 



(8) 



As it was discussed in Ref. [13[, most low-energy processes imply weak bounds on the coefficients 
xfj, which turn out to be of 0(1). However, some important constraints on tan/3 have started to 
appear, based on i?-physics |30t |. In order to discuss these results we find convenient to generalize 
the notation of Ref. [3V\ and define the couplings UidjH~^ and UidjH^ in terms of the matrices 



Xij, Yij and Zij (for leptons). In our case these matrices are given by: 

Xij = tan P 6ij 



sec P [rnd~ _ , 



Yii 



cot (3 5ii 



V2 

CSC 13 



m. 



~X.il 



where Xu and Yu are related with Sij and Pij defined in the Eq. ([8]) as follows: 



(9) 



Sij = ^ [(VcKM)jZ f^di Xij + rriu^ Yii{VcKM)lj\ 
1=1 
3 

Pij = ^ [(VcKM)i« rridi Xij - niu, YiiiVcKM)ij] 
1=1 



(10) 



The 33 elements of these matrices reduce to the expressions for the parameters X,Y,Z ( 



X33, 133, Z33) used in Ref.[31]. Based on the analysis oi B ^ Xg^ |3l|,l34], it is claimed that X < 20 
and Y < 1.7 for mjj+ > 250 GeV, while for a lighter charged Higgs boson mass, mjj+ ~ 200 GeV, 
one gets {X,Y) < (18,0.5). Fig. [1] shows the values of {X,Y) as a function of tan/3 within our 
model. Thus, we find important bounds: Ixssl ~ 1 for 0.1 < tan/3 < 70. Although in our model 
there are additional contributions (for instance from c-quarks, which are proportional to X23), 
they are not relevant because the Wilson coefficients in the analysis oi B ^ Xg'j are functions 



of m'^/M^r or m^/m|^+ [33l], that is, negligible when compared to the leading X33 effects, whose 



Wilson coefficients depend on m^/M^ or m1/m?^^. Other constraints on the charged Higgs mass 
and tan /3, based on Aa^j, the p parameter, as well as B-decays into the tau lepton, can be obtained 



34], 



351 ]. For instance, as can be read from Ref. [361] . one has that the decay B tu, implies a 
constraint such that for raij+ = 200 (300) GeV, values of tan/3 less than about 30 (50) are still 
allowed, within MSSM or 2HDM-III: However, these constraints can only be taken as estimates, 
as it is likely that they would be modified for 2HDM-III. In summary, we find that low energy 
constraints still allow to have xfj = 0{1)^ . 



A more detailed analysis that includes the most recent data is underway [37 
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FIG. 1: The figure shows X33, 133, X23 and F23 vs. tan/?, taking Xs's = 1 (solid), Xsli =0.1 (dashes), 
= -0-1 (dots), xs.a = -1 (dashes-dots). 



III. DECAYS OF THE CHARGED HIGGS BOSON 



Let us now discuss the decay modes of the charged Higgs boson within our model. Hereafter, 
we shall refer to four benchmark scenarios, namely, (i) Scenario A: xfj = 1, xfj = 1; (ii) 
Scenario B: xfj = O-l; xfj = 1; (iii) Scenario C: xfj = !> xfj = 0-1; (iv) Scenario D: x^j = 0.1, 

xfj = 0.1. We have performed the numerical analysis of charged Higgs boson decays by taking 
tan/3 = 0.1, 1, 15, 70 and varying the charged Higgs boson mass within the interval 100 GeV 
< mfj± < 1000 GeV, further fixing nij^o = 120 GeV, m^o = 300 GeV and the mixing angle at 
a = 7r/2. 

The condition ^-Si < i in the frame of the 2HDM-II implies ^"'^ ^ SGprn^ ^j^igj^ leads 

to 0.3 < tan/3 < 130. However, in the 2HDM-HI we have that ^ « t^^^ ( t^) > we 

^ cos a ^ 

have checked numerically that this leads to 0.08 < tan/3 < 200 when 1x^3 1 « 1 and 0.3 < tan/3 < 



130 as long as IX33I ~^ recovering the result for the case of the 2HDM-II [7], |38[]. In this sense, 
if we consider the constraints imposed by the perturbativity bound, a portion of the low tan/3 
appearing in some graphs would be excluded. However, we have decided to keep that range both 
to show the behaviour of the quantities of interest, and also because we have to keep in mind that 
such criteria (perturbativity) should be taken as an order of magnitude constraint. 

The expressions for the charged Higgs boson decay widths Uidj are of the form: 

T(H'^ ^ Uidj) = ^ — o — A^/^(m?r+,m,^,.,mj.) 



1 

X 

2 



2 22 



(Sf, + 4) - mu^m,^ (Sf^ - 4) ) , (11) 



where A is the usual kinematic factor A(a, b, c) = {a — b — c)^ — 46c. When we replace Xud — > 0, the 
formulae of the decays width become those of the 2HDM-II: see, e.g., Ref. Qj. Furthermore, the 
expressions for the charged Higgs boson decay widths of the bosonic modes remain the same as in 
the 2HDM-II. Then the results for the Branching Ratios (BRs) are shown in Figs. 2-8, and have 
the following characteristics. 

Scenario A. In Fig. we present the BRs for the channels — > tb, cb, ts, t^Vt-, W^hP , 
^ as a function of mu+, for tan/3 = 0.1 and fixing m/jO = 120 GeV, m^o = 300 GeV and the 
mixing angle a = 7r/2. When mij+ < 175 GeV , we can see that the dominant decay of the charged 
Higgs boson is via the mode c5, with BR(ffj^ ^ cb) ~ 1, which will have important consequences 
for charged Higgs boson production through c6-fusion at the LHC and may serve as a distinctive 
feature of this model. For the case 175 GeV < mjj+ < 180 GeV the mode ts is relevant, which it 
is also very different from the 2HDM-II and becomes an interesting phenomenological consequence 
of the 2HDM-III. We can also observe that, for mjj+ > 180 GeV, the decay mode tb is dominant 
(as in the 2HDM-II). Now, from Fig. [2]^b), where tan/3 = 1, we find that the dominant decay mode 
is into T^Ur for the range < 175 GeV, again for 175 GeV < mfj+ < 180 GeV the mode ts is 
the leading one, but for 180 GeV< mH+ < 600 GeV, the decay channel W+h° becomes relevant, 
whereas for the range 600 GeV < mfj+ the mode VK^A'^ is dominant. It is convenient to mention 
that this sub-scenario is special for the mode tb, because its decay width is zero at the tree-level, 
since the CKM contribution is canceled exactly with the terms of the four-zero texture implemented 
for the Yukawa coupling of the 2HDM-III. Then, see Fig. [2)[|c), for the case with tan/3 = 15 one gets 
that BR{H+ T+Vr) ^ 1 when 712^+ < 180 GeV. However for 180 GeV< 171^+ < 300 GeV, the 
dominant decay of the charged Higgs boson is the mode tb, while in the range 300 GeV< mff+ , the 
decay channel W~^h^ is also relevant. For tan/? = 70, we show in plot Fig. [2][|d) that the dominant 
decay of the charged Higgs boson is the mode t'^i't, when mj{+ < 300 GeV, but that, for 300 
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FIG. 2: The figure shows the BRs of the H'^ decaying into the principal modes in Scenario A, taking 
Xf. = 1, xf. = 1, m^o = 120 GeV, m^n = 300 GeV and a = tt/2 for: (a) tan/? = 0.1, (b) tan/3 1, (c) 
tan/3 = 15, (d) tan/3 = 70. The hnes in each graph correspond to: (1) BR{H+ tb), (2) BR{H+ cb), 
(3) BR{H+ ts), (4) BR{H+ r+i^r), (5) BR{H+ -> T^+Zi"), (6) BR(i?+ ^ W+A°). 

GeV< mjj+ < 400 GeV, the decay channel tb becomes the leading one, whereas for the range 400 
GeV < mfj+, the mode W~^h^ is again dominant. 

Scenario B. In Fig. [3l we present the BRs of the channels — > tb, cb, ts, t'^i't, W'^h'^, W~^A^ 
as a function of mjj+. From Fig. E^a), we observe that for tan (3 = 0.1, when mfj+ < 175 GeV, the 
dominant decay of the charged Higgs boson is the mode cb, with BIi{H^ cb) ~ 1. When 175 
GeV < mjj+ < 180 GeV the mode ts is important and for mjj+ > 180 GeV the decay mode tb 
becomes the leading one. From Fig. El^b), we see that, for tan/3 = 1, the dominant decay mode is 
now into r^i'r for mjj+ < 175 GeV, while in the range 175 GeV < mi[+ < 180 GeV the mode ts is 
relevant. For 180 GeV< mjj+ < 500 GeV the decay channel tb becomes the leading one, whereas 
for the range 500 GeV < mfj+ the mode VF^^'' is dominant. From Fig. [3|^c), with tan/3 = 15, one 
gets that BK{H^ — > t^i^t) ~ 1 for mi{+ < 180 GeV. For 180 GeV< m^+, the dominant decay of 
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FIG. 3: The same as in Fig.[2]but taking Xy- = 0.1, xfj = 1 (Scenario B). 

the charged Higgs boson is instead the mode W~^h^. Then, for tan/3 = 70, we show in Fig. [3|^d) 
that the dominant decay of the charged Higgs state is via the mode t^Ut when mj^+ < 350 GeV 
while for 350 GeV< mf{+ the decay channel W~^h^ becomes the leading one. 
Scenario C. In Fig. [5] we show the corresponding plots for the BRs of the channels — > tb, cb, 
ts, T^fri W'^hP^ as a function of mfj+. For tan /3 = 0.1, as one can see in Fig. [D^a), the 

mode cb is dominant when mff+ < 170 GeV, but for 175 GeV < mjj+ < 180 GeV the mode ts 
is relevant, while for 180 GeV< mff+ the mode tb becomes dominant. For tan/3 = 1, we observe 
from Fig. IHb) that the dominant decay modes are: t'^i't in the range mjj+ < 170 GeV, ts for 
175 GeV < mH+ < 180 GeV, W+h^ for 180 GeV < mH+ < 600 GeV and when 600 GeV 

< mff+. For tan/3 = 15, as shown in Fig. IDj^c), the relevant decay channels are: t^i't in the 
range < 180 GeV, tb when 180 GeV < mH+ < 300 GeV, for 300 GeV < 771^+. In 

Fig. Hl^d), for tan/3 = 70, we observe that t^Ut dominates when mjj+ < 180 GeV, but when 180 
GeV < 171^+ < 900 GeV the mode tb is the leading one, whereas for 900 GeV < mij+ the mode 
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FIG. 4: The same as in Fig.H but taking = 1, xfj = 0.1 (Scenario C). 
W~^h^ is the most relevant one. 

Scenario D. In Fig. [5] we present plots for the BRs of the channels tb, cb, ts, t~^Vt, W~^h^, W~^A^ 
as a function of mjj+. For tan/3 = 0.1, we show in Fig. [S][^a) that the dominant decay modes for 
the are: cb in the range 111^+ < 175 GeV, ts when 175 GeV < 111^+ < 180 GeV, tb for 180 
GeV < mjj+. For tan/3 = 1, we show in Fig. [5||b) that the mode r'^i'r is dominant in the range 
mfj+ < 175 GeV, whereas for 175 GeV < mjj+ < 180 GeV the relevant decay channel is ts, whilst 
the mode tb dominates for 180 GeV < mf{+ < 550 GeV and the mode W^^^*^ does so when 550 
GeV < mjj+. For tan/3 = 15, we observe in Fig.[5|^c) that the relevant decay channels are: t'^Vt in 
the range mH+ < 250 GeV and W~^h° for 250 GeV < mH+- Finally, for tan/3 = 70, see Fig. [5{d), 
we obtain that, when mjj+ < 230 GeV, the mode t'^i't becomes the most important one but, for 
230 GeV < mfj+ < 800 GeV, the channel tb is the leading one, whereas, for 800 GeV < mfj+, the 
mode W~^h^ is the dominant one. 
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FIG. 5: The same as in Fig. El but taking x", = 0.1, xi = 0.1 (Scenario D). 



In order to cover further the Higgs sector in our analysis, it is appropriate to also mention 
how the previous results change with m/jO, m^o and a. Regarding the former two, clearly, the 
lather the neutral Higgs boson mass the later the corresponding decay channel will onset. 
Regarding the latter, we adopted two further choices, a = (3 and 0, in all scenarios previously 
studied. In general the behavior of the decay modes of the charged Higgs boson is similar to the 
cases presented above, except for the decay channel Wh^. For a = 0, this mode has BR < 10^'^ 
when tan/3 is large. However, for tan /? < 1, it becomes the dominant one. In the case a = f3, the 
decay channel Wh^ can be dominant with a BR that could be 0{1). 

As a general lesson from this section, and distinctive features of our 2HDM-III, we can see 
that both decay modes W^h^ and cb become very relevant phenomenologically, effectively of 0(1) 
for some of the scenarios considered. Therefore, we want to study next the general behaviour of 
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FIG. 6: The figure shows the BR{H+ c6)iii/BR(i7+ cb)u vs. mH+, taking tan/? = 0.1, 1, 15, 70 for 
{^j(:f = 1, (b) x:tf = -1, (c)x:tf = 0.1, (d) = -0.1. 

these decay modes, in relation to the 2HDM-II case. In order to compare the 2HDM-III results 
with those in the 2HDM-II, we show in Fig. [6] the ratio BR(F+ c6)ni/BR(i7+ ch)ii vs. 
mH+, taking again tan/3 = 0.1, 1, 15, 70, for: (a)x"/ = 1, (b) x^/ = -1> (c) X^f = 0.1 and 
(d) Xij'^ = —0.1. We observe that the mode cb is important when 200 GeV < mjf+ < 300 GeV 
and for 0.1 < tan/3 < 1, taking Xij'^ = 1- Now, in Fig. [71 we present the behaviour of the ratio 
BR(i7+ — > cb)iu/BK{H^ cb)u as a function of xtf: for tbe cases: (a) mjf+ = 150 GeV, (b) 
mfj+ = 300 GeV, (c) mfj+ = 450 GeV and (d) mfj+ = 600 GeV. Again, one can see that the 
largest enhancement arises when 'mj^+ = 300 GeV and x^j'^ = 1- Finally, specific to the 2HDM-III, 
we show in Fig. [8] the ratio BR(iJ+ W+h%u/BR{H+ cb)ui vs. mH+, taking tan/3 = 0.1, 
1, 15, 70, for: (a)^^/ = 0.1, (b) x^/ = "d, (c) X^i" = 1 and (d) x-f = 1- We find that 
BR(i7^ l^^^*^):!! is much larger than BR(i7^ cb)iu when Xij'^ = 1 and the mass of the 
charged Higgs boson is close to the upper limit obtained by unitarity conditions, which is about 
800 GeV [2u\. Thus, we find that the effect of the modified Higgs couplings typical of the 2HDM-III 
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FIG. 7: The figure shows the BR{H+ c6)iii/BR(iJ+ -> cb)n vs. Xif, taking tan/3 = 0.1, 1, 15, 70 for 
(a)TOH+ = 150 GeV, (b) mH+ ^ 300 GeV, (c)mff+ = 450 GeV, (d) tojj+ = 600 GeV. 

shows up clearly in the pattern of charged Higgs boson decays, which can be very different from 
the 2HDM-II case and thus enrich the possibilities to search for states at current (Tevatron) 
and future (LHC, ILC/CLIC) machines. 



IV. CHARGED HIGGS BOSON PRODUCTION AT LHC 



The production of charged Higgs bosons at hadron colliders has been evaluated in early 



(also for the Superconducting Super Collider, SSC) and more recent [22| (for the LHC) literature, 
mainly for the 2HDM-II and its SUSY realization (i.e., the MSSM). In these two scenarios, when 
kinematically allowed, the top quark decay channel t bH^ is the dominant production 
mechanism. Instead, above the threshold for such a decay, the dominant production reaction 
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FIG. 8: The figure shows the BR{H+ W+h°)ni/BR{H+ cb)m vs. mH+, taking tan/3 = 0.1, 1, 15, 70 
for {a)xlj'^ = 0.1, (b) x^/ = -O-l. Wx^/ = 1> (d) X^/ = "l- 

is gluon-gluon fusion into a 3-body final state, i.e., gg —>■ tbH^^. Both processes depend on the 
coupHng H~tb and are therefore sensitive to the modifications that arise in the 2HDM-III for this 
vertex. However, detection of the final state will depend on the charged Higgs boson decay mode, 
which could include a complicated final state, that could in turn be difficult to reconstruct. For 
these reasons, it is very important to look for other production channels, which may be easier 
to reconstruct. In this regard, the s-channel production of charged Higgs bosons, through the 
mechanism of c6-fusion, could help to make more viable the detection of several charged Higgs 
boson decay channels jg?!. 

Here, we shall evaluate the predictions of the 2HDM-III for the t — > hH^ (and sH^) decay rate 
plus the cb- as well as the gg-fusion mechanisms (hereafter, referred to as 'direct' and 'indirect' 
production, respectively). 



In fact, these two mechanisms are intimately related, see below. 
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A. The decays t ^ H+ b, H+s 

We shall discuss here the charged Higgs boson interactions with heavy quarks (t, b, c, s) and their 
implications for charged Higgs boson production through top quark decays. In order to study the 
top quark BRs, besides the SM decay mode t — > bW~^, we need to consider both decays t — > bH^ 
and t — > sH^, because these modes could both be important for several parameter configurations 
within our model. The decay width of these modes takes the following form: 

n,2 



lioTrm^m^ 



H+ 



SI, + 



(mt + TTife)^ — m 



H+ 



Pi , (12) 



where A is the usual kinematic factor A(a, b,c) = {a — b — cf' — Abe, j = 2 for the mode sH^ and 
j = 3 for the mode bH^ . Furthermore, we shall neglect the decay width for the light fermion 
generations. If one takes XiJ ~^ 0, the formulae for the decay width reduce to the 2HDM-II case: 
see, e.g., ^J. 

We have explored several theoretically allowed regions within our scenario, which are constrained 



by using the bounds on the BR(t bH^). In the so-called "tauonic Higgs model" 23|, the decay 



mode {H^ u-j-) dominates the charged Higgs boson decay width, and BR(t bH^) is 



constrained to be less than 0.4 at 95 % C.L. 



23l |. However, if no assumption is made on the 



charged Higgs boson decay, BR(t bH^) is constrained to be less than 0.91 at 95 % C.L. \23l\ . 
However, the combined LEP data exclude a charged Higgs boson with mass less than 79.3 GeV 
at 95 % C.L., a limit valid for an arbitrary BR(/7+ I'r) [25]. Thus, in order to perform our 

analysis, we need to discuss all the charged Higgs boson decays following the steps of our previous 
paper 



19l |. In the present section, we take all charged Higgs boson decays relevant for masses below 
that of the top quark, thus including the modes T^iyr,ts,cb,W~^h^ ,W~^A^. As usual, we refer to 
our four benchmark scenarios. 

Scenario A. Remember that this scenario was defined by taking xfj = 1 and xfj = while for 
tan/? we considered the values tan /? = 0.1, 1, 15, 70. In Fig. [9] we present plots of BR(t ^ 6 us. 
m'fj and BR(t s H^) vs. mfj+ . We can observe that a charged Higgs boson within the mass range 
80 GeV < mH+ < 170 GeV and for 1 < tan/3 < 70 satisfies the constraint BR{t bH+) < 0.4. 
Furthermore, from the plots of Fig. [21 we can see that in this scenario the dominant decay mode is 
into t'^Vt for tan /? = 1, 15, 70, therefore we fall within the scope of the tauonic Higgs model, so that 
BR(t — > b) < 0.4 applies. However, for the case tan /? = 0.1, the dominant decay of the charged 
Higgs boson is the mode cb and the mode t — > bH'^ satisfies the constraint BR(t — > bH^) < 0.9 
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FIG. 9: It is plotted: the BR{t bH+) vs. r7i^+ (left); the BR(t s H+) vs. mH+ (right). Here is for 
Scenario A, obtained by taking Xij = 1 and xfj = li for tan/3 = 0.1 (solid), 1 (dashes), 15 (dots) and 70 
(dashes-dots). 

in the range described above. 

Scenario B. In Fig. [TO]we present similar plots for the case Xij = 0-1 ^iid xfj = 1; taking tan j3 = 
0.1, 1, 15, 70. We can observe that the mode t — > bH~^ satisfies the constraint BR(t hH^) < 0.4 
within the ranges 80 GeV < mjj+ < 170 GeV and 1 < tan/3 < 70. Thus, from Fig. El we can see 
that in this range the dominant decay mode is into r^Ur^ therefore this setups also falls within the 
realm of the tauonic Higgs model, so that BR(t h) < 0.4 must hold in this scenario. For 

tan P = 0.1, the dominant decay of the charged Higgs boson is cb, thus the channel t — > bH^ must 
satisfy the constraint BR(t bH^) < 0.9, which is fulfilled in the range studied. 
Scenario C. In Fig. [11] we present the corresponding plots for the case xfj = 1 xfj = 0-1) 
taking again tan /3 = 0.1, 1, 15, 70. We can observe that the mode t bH^ satisfies the constraint 
BR(t bH+) < 0.4 in the range 80 GeV < mj^+ < 170 GeV and 1 < tan/3 < 70. Similarly, as 
in scenario A, for tan/3 = 0.1 the dominant decay of the charged Higgs boson is the mode c6, thus 
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FIG. 10: The same as in Fig.lHlbut taking Xy- = 0.1 and xfj = 1 (Scenario B). 

the mode t bH^ must satisfies the constraint BR(t — > hH^) < 0.9, indeed satisfied in the range 
analyzed here. 

Scenario D. Recah that this was defined by taking x'^j = 0.1 and xfj = O-l- Iii Fig. [T2] we 
present the usual plots of the BR(t bH^) and BR(t bH'^) vs. mjj+. One can see that, for 
charged Higgs boson masses within the range 80 GeV < < 170 GeV and 1 < tan/3 < 70, 

the model fulfills the constraint BR(t — > bH^) < 0.4. Furthermore, for tan/3 = 0.1, the dominant 
decay of the charged Higgs boson is the mode cb, thus the mode t bH^ satisfies the constraint 
BR(t bH^) < 0.9 in the range studied. 

In short, as bottomline of these execises, we have identified regions of the 2HDM-III parameter 
space where a charged Higgs mass below mt — mi, has not been excluded by Tevatron. Therefore, 
the LHC is best positioned in order to probe charged Higgs bosons with such masses. 

B. Direct production of charged Higgs bosons at the LHC 



The H^qq' vertex with large flavor mixing coupling, that arises in the 2HDM-III, enables 
the possibility of studying the production of charged Higgs boson via the s-channel production 
mechanism, cb H'^ + c.c. This process was discussed first by Ref. 27|], both within topcolor 
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FIG. 11: The same as in Fig. [Slbut taking x", — 1 and xfj = 0-1 (Scenario C). 
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FIG. 12: The same as in Fig.[9]but taking xf, = 0.1 and xi = 0.1 (Scenario D). 
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FIG. 13: The figure shows the total cross section rates of process hih2{cb) H^X as a function of m//+ 
in the 2HDM-III at LHC energies (^/^ = 14 TcV), by taking xfs, = 1 and = 1 = 1,2,3). The lines 
correspond to: tan/3 = 0.1, tan/3 = 1, tan/3 = 15, tan/3 = 70. 

models and a simplified version of the 2HDM-III. Then the SUSY case was discussed in [sO^ and j4o| . 
Here we perform a detailed study of this mechanism within the 2HDM-III, paying special attention 
to the effects induced by the assumed Yukawa texture on the charged Higgs boson couplings. 
Defining the H^qq' coupling here as Cl^—^ + C/j^-y^, we can express the total cross section for 
direct production at hadron colliders as 



a{h,h2{d>) - H+X) = + /'f ^ (13) 

where 

' J r 



and r = mj^±/s. The Parton Distribution Functions (PDFs) fq^ix^Q^) used here are from 41], 
with scale choice = . 

From Eq. ([6]) we see that, for the case of the 2HDM-III, Cl and Cr entering the subprocess 
cb are given by 
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FIG. 14: The figure shows the total cross section rates of process hih2{cb) H^X as a function of mjLf+ 
in the 2HDM-II at LHC energies (^s = 14 TeV), by taking (Vckm)23 = 4.16 x 10"^ and (VbKM)33 ~ 1- 
The hnes correspond to: tan/3 = 0.1, tan/3 — 1, tan/3 — 15, tan/3 = 70. 



and 



V2Mw iri 



'o X CSC/3 

cot 13 rUc 021 - ^/rricmu, X21 



(Vckm)«3) 



Cr^C 



III 

R 



^ J. sec/3 , J 



CKMj2«- 



(15) 



(16) 



We notice here that Eqs. ([15]) and (fT6|) reduce to the case of the 2HDM-II if one takes 



Cl ^Cl = - -jJ— cot /? (ycKM)23 (17) 



and 



= &k = -^i^ tan/3mb(ycKM)23. (18) 
V 2Mw 

In Fig. [T3l we present plots for the total cross section rates of process hih2{cb) H^X as a 

function of mff+ in the framework of the 2HDM-III, by taking xfs = 1 s-i^cl X21 = 1 (^ = 1;2,3), 
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at LHC energies {^/s = 14 TeV), for the cases: (a) tan/3 = 0.1, (b) tan/3 = 1, (c) tan/3 = 15, 
(d) tan f3 = 70. The sum over I is performed over ah the three quark famiUes and we take for the 
quark masses: m„ = 2.55 MeV, = 5.04 MeV, nic = 1.27 GeV, = 104 MeV, ruh = 4.20 
GeV, and mt = 171.2 GeV [25]. We have checked numerically that the term proportional to 
^ csc^ (3 rricmt Ixfs (^CKM)33p provides the most important contribution to the cross section rates 
and dominates by far for X23 ^ 1- the other hand, the expected integrated luminosity at LHC 
is of the order 10^ pb~^ and given that a > 10~^ pb even for tan/3 = 70 and < 600 GeV, we 
can conclude that, in the context of the 2HDM-III, it is likely that a charged Higgs boson could 
be observed at LHC energies by exploiting direct production. 

In Fig. [T^ we present results for the total cross section rates of process /ii/i2(c6) H^X as a 
function of in the 2HDM-II at LHC energies {y/s = 14 TeV), by taking (V'ckm)23 = 4.16 x 10"^ 
and (Vckm)33 ~ 1, for the cases: (a) tan/3 = 0.1, (b) tan/3 = 1, (c) tan/3 = 15, and (d) tan/3 = 70. 
As we have already said, the expected integrated luminosity at LHC is of the order 10^ pb~^, hence 
we also conclude from this figure that in the framework of the 2HDM-II we obtain production rates 
for the charged Higgs boson via c6- fusion that may be detectable at LHC energies. 

C. Indirect production of charged Higgs bosons at the LHC 

We have found that, in some of the 2HDM-III scenarios envisaged here, light charged Higgs 
bosons could exist that have not been excluded by current experimental bounds, chiefly from LEP2 
and Tevatron. Their discovery potential should therefore be studied in view of the upcoming LHC 
and we shall then turn our attention now to presenting the corresponding hadro-production cross 
sections via an indirect channel, i.e., other than as secondary products in (anti)top quark decays 
and via c6-fusion, considered previously. 

As dealt with so far, if the charged Higgs boson mass mij± satisfies mjj± < mt — rrib, where 
mt is the top quark mass and m^ the bottom quark mass, particles could be produced in the 
decay of on-shell (i.e., Tt 0) top (anti-)quarks t — > bH^ and the c.c. process, the latter being in 
turn produced in pairs via qq annihilation and gg fusion. We denote such a production channel 
as qq, gg ^ tt ^ tbH^ + c.c. (i.e., if due to (anti-)top decays) whilst we use the notation qq, 
gg thH^ + c.c. to signify when further production diagrams are included^. In fact, owing to the 

^ Altogether, they represent the full gauge invariant set of Feynman graphs pertaining to the 2 — > 3 body process 
with a thH~ + c.c. final state: two for the case of qq annihilation and eight for gluon-gluon fusion, see, e.g., 
Eq. (1.1) of Hi]. 



23 




FIG. 15: The figure shows the cross sections of production at the LHC through the channel qq,gg 
tbH^ + c.c. in Scenario A {xfj = 1 and xfj = 1) and for tan/? = 0.1, 1, 15 , 70. 



large top decay width (Tt > 1.5 GeV) and due to the additional diagrams which do not proceed 



43, |4J, l45| , charged Higgs 



via direct tt production but yield the same final state tbH~ + c.c. 
bosons could also be produced at and beyond the kinematic top decay threshold. The importance 
of these effects in the so-called 'threshold' or 'transition' region (mu± ~ mt) was emphasized in 



various Les E 
of Refs. 



HERWIG 



48l . l49l . 150( 1 , so that the calculations 



441 (based on the appropriate g^gg — > tbH^ description) are now implemented in 



ouches proceedings [46|, |47l] as well as in Refs 
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and PYTHIA [551, 



5g]. A comparison between the two generators was 



carried out in Ref. [481]. For any realistic simulation of production with mfj± ^ mt, as can well 
be the case here, the use of either of these two implementations is of paramount importance. 

Here, we use HERWIG version 6.510 in default configuration, by onsetting the subprocess 
IPROC = 3839, wherein we have overwritten the default MSSM/2HDM couplings and masses with 
those pertaining to the 2HDM-ni: see Eqs. ©-([H]). The production cross sections are found in 
Figs. [TSHlHl for our usual scenarios: A {xfj = 1 and xfj = l)i B {xfj = 1 and xfj = 0-1)) C 



{xfj = 0.1 and xfj = 1) a-iid D {xfj = 0.1 and xfj = 0.1) As usual, we adopt our four choices of 



tan (3. 
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FIG. 17: The same as in Fig. [15] but taking x^- = 0.1 and xf, = 1 (Scenario C). 
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FIG. 18: The same as in Fig.[T5]but taking x% = 0.1 and xtj = 0.1 (Scenario D). 

Altogether, by comparing the q q, g o — > tbH~ + c.c. cross sections herein with, e.g., those of the 
MSSM in [63 or the 2HDM in [49|, l57[, it is clear that the 2HDM-III rates can be very large and 
thus the discovery potential in ATLAS and CMS can be substantial, particularly for a very light 
H^, which may pertain to our 2HDM-III but not the MSSM or 2HDM-II. However, it is only by 
combining the production rates of this section with the decay ones of the previous ones that actual 
event numbers at the LHC can be predicted. 



V. EVENT RATES OF CHARGED HIGGS BOSONS AT THE LHC 

To illustrate the type of charged Higgs signatures that have the potential to be detectable at the 
LHC in the 2HDM-ni, we show in Tabs. HI and HIl the event rates of charged Higgs boson through 
the channels qq,gg — > tbH^ + c.c. and cb — > + c.c, alongside the corresponding production 
cross sections (c's) and relevant BRs, for a combination of masses, tan /? and specific 2HDM-ni 
parameters amongst those used in the previous sections (assuming m/jO = 120 GeV, m^o = 300 
GeV and the mixing angle at a = tt/2 throughout). In particular, we focus on those cases where 
the charged Higgs boson mass is above the threshold for t bH'^ , for two reasons. On the one 
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hand, the scope of the LHC in accessing t — > bH~^ decays has been estabhshed in a rather model 
independent way. On the other hand, we have dealt at length with the corresponding BRs in 
section III. (As default, we also assume an integrated luminosity of 10^ pb^'^.) 

To illustrate these results, let us comment on one case within each scenario. From TablelH we can 
see that for Scenario A, with (x^- = 1, xfj = 1) and tan/3 = 15, we have that the is heavier than 
mt — mfj, as we take a mass mf^+ = 400, thus precluding top decay contributions, so that in this case 
a{pp tbH^) 2.2 x 10~^ pb, while the dominant decays are tb, t'^v^ W^h!^, W^A^ which 

give a number of events of 7040, 46, 13860, 374, respectively. In this case the most promising signal 
is W^h?. However, when tan/? = 70 we have that all event rates increase substantially. 

Here, the signal W'^h^ is still the most important with an event rate of 15480. 

Then, for Scenario B (x^^. = 0.1, xfj = 1), we have that is again above the threshold 
for t H^b. So, for the declared values of the relevant parameters, we take a charged Higgs 
boson mass mjj+ = 600 for tan/3 = 1 and tan/3 = 70, respectively. In such a case the decay 
W^h^ can reach significant numbers for the LHC. We obtain a number of events of 3960 
and 2703, respectively. The other decay that has a large BR is — > M/^"*"j4'', and in these cases 
the number of events ranges over 1200-3500. 

Next, we discuss the Scenario C {xfj = ^^xfj = 0-1) for tan/3 = 15. Here, we obtain that the 
signals H'^ tb and W~^h^ are the most relevant ones, with a number of events about 34560 and 
26240, respectively. 

Finally, for scenario D {xfj = 0.11, xfj = 0.1) the dominant decays are — > tb,T^i>r and 
W~^h^, which give a spectacular number of events: 269800, 68400 and 34200, respectively. Here, 
we have set tan /3 = 70. 

All these rates correspond to the case of indirect production. The contribution due to direct 
production is in fact subleading, especially at large mjj± values. Nonetheless, in some benchmark 
cases, they could represent a sizable addition to the signal event rates. This is especially the case 
for Scenario A with tan /? = 15 or 70 and Scenario C with tan /? = 15. In general though, also 
considering the absence of an accompanying trigger alongside the H^,i.e. for instance a top quark 
produced in gb H^t could help to identify the signal. Thus, we expect that the impact of 
c5-fusion at the LHC will be more marginal that that of gg-lusion for large Higgs masses, in fact, 
at times even smaller that the contribution from gg-annihilation. 
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TABLE I: Summary of LHC event rates for some parameter combinations within Scenarios A, B, C, D with 
for an integrated luminosity of 10^ pb~^, for several different signatures, through the channel qq, gg — > tbH^ 
+ c.c. 





tan fS 


Tnf{+ in GeV 


a{pp — » H'^tb) in pb 


Relevant BRs 


Nr. Events 


(1,1) 


15 


400 


2.23 X 10-1 


BR(_H'+ ^ tb) « 3.2 X 10-1 
BR(_H'+ T+u^) « 2.1 X 10-3 
BR(_H'+ W+h°) « 6.3 X lO-i 
BR(_ff+ W+A"^ « 1.7 X 10-2 


7040 
46 
13860 
374 


(1,1) 


70 


400 


4.3 X 10-1 


BR(if+ tb) « 3.5 X 10-1 
BR{H+ cb) Si 1.4 X 10-2 

BR(_ff+ T+U-r) w 2.5 X 10"! 

BR(_ff+ VK+ftO) Ki 3.6 X 10-1 


15050 

602 
10750 
15480 


(0.1,1) 


1 


600 


1.1 X 10-1 


BR(_f/+ ^ tb) Si 3 X 10-1 
BR(if+ -* ts) « 9.1 X 10-* 

BR(_ff+ W+h'^) ^ 3.6 X 10-1 
BR(_ff+ -» W+AO) ^ 3.2 X 10-1 


3300 
10 

3960 
3520 


(0.1,1) 


70 


600 


5.1 X 10-2 


BR(_ff+ T+Ur) « 1.2 X 10-1 
BR(_ff+ ^ tb) « 9.4 X 10-2 
BK{H+ W+hP) » 5.3 X lO-i 
BK{H+ W+A°) Ri 2.3 X lO-i 


612 
470 
2703 
1173 


(1,0.1) 


15 


300 


6.4 X 10-1 


BR(_ff+ -> tb) Ri 5.4 X 10-1 
BR(_ff+ ^ c6) Si 5 X 10-1 
BR(_ff+ T+Ur) 3.9 X 10-2 

BK(yi 1 — ir 1 /)») % i.i X 10 1 


34560 
32 
2535 

2()2 10 


(0.1,0.1) 


70 


300 


3.8 


BR(//+ T+Vr) ^ 1.8 X IQ-i 
BR(//+ tb) ^ 7.1 X 10-1 
BR(//+ cb) « 2.4 X 10-3 
BR(H+ -» W+h°) « 9 X 10-2 


68400 
269800 
912 
34200 



VI. CONCLUSIONS 



We have discussed the implications of assuming a four-zero Yukawa texture for the properties 
of the charged Higgs boson, within the context of a 2HDM-in. In particular, we have presented 
a detailed discussion of the charged Higgs boson couplings to heavy fermions and the resulting 
pattern for its decays. The latter clearly reflect the different coupling structure of the 2HDM-III, 
e.g., with respect to the 2HDM-II, so that one has at disposal more possibilities to search for 
states at current and future colliders, ideally enabling one to distinguish between different Higgs 
models of EWSB. We have then concentrated our analysis to the case of the LHC and showed 
that the production rates of charged Higgs bosons at the LHC is sensitive to the modifications 
of the Higgs boson couplings. We have done so by evaluating 2HDM-HI effects on the top decay 
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TABLE II: Summary of LHC event rates for some parameter combinations within Scenarios A, B, C, D with 
for an integrated luminosity of 10^ pb~^, for several different signatures, through the channel c6 — > + 
c.c. 





tan /3 


mu-\- in GeV 


a{pp H+ + X) in pb 


Relevant BRs 


Nr. Events 










BR(H+ tb) Ri 3.2 X 10-1 


3648 


(1,1) 


15 


400 


1.14 X 10-1 


BR{H+ T+u°) ^ 2.1 X 10--^ 
BR{H+ W+h") ^ 6.3 X lO-i 
Br(^H+ W+A^^ « 1.7 X 10-2 


24 
7182 
194 










BR(i?+ tb) « 3.5 X 10-1 


4375 


(1,1) 


70 


400 


1.25 X 10-1 


BR{H+ ch) ^ 1.4 X 10-2 
BR(/i'+ T+u^) Ri 2.5 X 10-1 
BR(//+ W+h'^) ^ 3.6 X 10-1 


175 
3125 
4500 










BR(//+ ^tb)^3x 10-1 


10 


(0.1,1) 


1 


600 


3.41 X 10-4 


BR(i?+ ts) « 9.1 X lO--* 
BR{H+ W+hO) X 3.6 x lO-i 
BR{H+ W+A°) « 3.2 X lO-i 



12 

11 










BR{H+ -> T+u^) « 1.2 X 10-1 


24 


(0.1,1) 


70 


600 


1.98 X 10-^ 


BR(//+ tb) ^ 9.4 X 10-2 
BR{H+ W+hP) » 5.3 X lO-i 
BR{H+ -> W+A°) « 2.3 X lO^i 


19 
105 

45 










BR(//+ tb) ^ 5.4 X 10-1 


21546 


(1,0.1) 


15 


300 


3.99 X 10-1 


BR{H+ cb) ^5x 10-'' 
BR(Jf"+ T+i^r) « 3.9 X 10-2 
BR(if + W+h°) « 4.1 X 10-1 


20 
1556 
16359 










BR(_ff+ ^ T+Ur) « 1.8 X 10-1 


6984 


(0.1,0.1) 


70 


300 


3.88 X 10-1 


BR(_ff+ tb) ^ 7.1 X 10-1 
BR(_ff+ cb) « 2.4 X 10-3 
BR(_H"+ W+h°) « 9 X 10-2 


27548 
93 
3492 



t bH'^ as well as in the s-channel production of through c6-fusion and the multibody final 
state induced by gg-iusion and ^'^-annihilation. Finally, we have determined the number of events 
for the most promising LHC signatures of a belonging to a 2HDM-III, for both cb + c.c. 

and qq ibH^ + c.c. scatterings (the latter affording larger rates than the former). Armed with 
these results, we are now in a position to carry out a detailed study of signal and background rates, 
in order to determine the precise detectability level of each signature. However, this is beyond the 
scope of present work and will be the subject of a future publication. 
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